Population III Star Formation and IMF 



Michael L. Norman 

Center for Astrophysics and Space Sciences, University of California, San Diego, La Jolla, California USA; 

mlnorman @ ucsd. edu 

Abstract. We review recent 3D cosmological hydrodynamic simulations of primordial star formation from cosmological 
initial conditions (Pop III. 1) and from initial conditions that have been altered by radiative feedback from stellar sources (Pop 
III. 2). We concentrate on simulations that resolve the formation of the gravitationally unstable cloud cores in mini-halos over 
the mass range 10^ < M/M,;, < 10^ and follow their evolution to densities of at least 10'*'cm^-' and length scales of < 10^^ 
pc such that accretion rates can be estimated. The advent of ensembles of such simulations exploring a variety of conditions 
permits us to assess the robustness of the standard model for Pop III. 1 star formation and investigate scatter in their formation 
redshifts and accretion rates, thereby providing much needed information about the Pop III IMF. The simulations confirm 
the prediction that Pop III.l stars were massive (~ 100 Mm), and form in isolation in primordial mini-halos. Simulations of 
Pop III. 2 star forming in relic HII regions suggest somewhat lower masses (~ 3OM0) which may help explain the chemical 
abundances of extremely metal poor stars. We note that no 3D simulation at present has achieved stellar density let alone 
followed the entire accretion history of the star in any scenario, and thus the IMF of Pop III stars remains poorly determined 
theoretically. 
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1. THE STORY UP TO 2003 

Dark matter mini-halos in the mass range 10^ — 10^ solar 
mass viriaHzing at high redshift are believed to collect 
enough primordial gas within them to host the forma- 
tion of the first generation of stars (Population III). Me- 
diated by molecular hydrogen cooling, gas in the cen- 
ters of these halos condenses and becomes unstable to 
gravitational fragmentation with a typical mass scale of 
a few hundred solar masses. Pioneering high resolution 
3D hydrodynamic cosmological simulations by Abel et 
al. (2002) and Bromm et al. (2002), and more recently 
by Yoshida et al. (2006) and O'Shea & Norman (2007a) 
have shown that gravitational fragmentation produces 
only one massive fragment per halo. If this gas does not 
fragment further as it approaches stellar density. Pop III 
stars with this typical mass would be produced. The fates 
and feedback effects of such stars has many interesting 
consequences for the early structure formation as dis- 
cussed elsewhere in this volume (see review by Ciardi). 

While this picture is now firmly established, based as 
it is on rather simple physics elucidated by the simula- 
tions, it is important to realize the remaining uncertain- 
ties. First, as of 2003 when the first reviews of primor- 
dial star formation were being written (Barkana & Loeb 
2001, Bromm & Larson 2004, Glover 2005), a rather 
small number of such simulations had been done, rais- 
ing the question of how ubiquitous this mode of star for- 
mation was in the early universe. Second, despite the 
large range of scales covered by the simulations, they 
were necessarily terminated well before stellar density 



was reached due to missing physics. The possibility of 
sub-fragmentation could not be ruled out, with conse- 
quent uncertainty on the mass scale of the first stars. 
Third, as emphasized by Barkana & Loeb (2004), for 
technical reasons the simulations were carried out in 
quite small cosmological volumes (< 1 Mpc comoving), 
the result being that the redshift of formation of the first 
stars found (z ^ 20) was underestimated. Larger volumes 
would contain rarer peaks in the density field, and these 
would presumably form Pop III stars earlier (White & 
Springel 2000). Could the higher densities and tempera- 
tures of this earlier epoch alter the mass scale of the truly 
first stars? And finally, there was a host of complicated 
feedback effects-radiative, kinetic, and chemical-which 
were not included in these first simulations which could 
"mess up" the simple picture by altering the initial con- 
ditions and the cooling properties of the collapsing cloud 
cores. 

The most concerning of these feedback effects was 
the build-up of an FUV background by emissions from 
the first stars that could photo-dissociate the hydrogen 
molecules which mediate Pop III star formation in the 
first place (Haiman, Rees & Loeb 1997, Ciardi, Fer- 
rara & Abel 2000, Haiman, Abel & Rees 2001). The 
first simulations examining this "negative feedback ef- 
fect" (Machacek, Bryan & Abel 2001, hereafter MBA; 
and Yoshida et al 2003, hereafter Y03) concluded that 
Pop III star formation would not be suppressed, but 
merely delayed. The concept of a critical halo mass ca- 
pable of forming a Pop III star was introduced by MBA. 
They found that the critical halo mass is an increas- 
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ing function of the FUV background mean intensity in 
the Lyman-Werner bands (11.2-13.6 eV) suggesting that 
Pop III star formation would become self-hmiting. Y03 
found that gas cooling becomes inefficient above J21 = 
0.01 due to the photo-destruction of H2. Here J21 is the 
mean intensity of the FUV background in units of 10^^' 
ergs/cm^/sec/ster/Hz. Using this information Y03 con- 
structed a semi-analytic model of the cosmic history of 
Pop 111 star formation. They found a rapidly rising pop- 
ulation of Pop 111 stars over the redshift interval 35 > z 
> 20 "the rise of Pop 111", which begins to become self- 
regulated at z ~ 25 when the FUV background attains 
721 = 0.01 . They did not consider the lower redshift evo- 
lution as this was out of the range of their simulations, 
but rather speculated that Pop II star formation domi- 
nated the cosmic star formation history below z ^ 20. 

It is important to point out that neither the MBA simu- 
lations nor the Y03 simulations were evolved below red- 
shifts of ~ 20 or had the spatial resolution to follow 
the chemo-thermal-hydrodynamics of cloud collapse on 
scales below ~ 0.1 pc, making "the fall of Pop III" by 
radiative feedbacks alone quite uncertain and understud- 
ied. Indeed, Y03 found the Pop III global SFR was still 
rising at the end of their simulations, implying that we 
don't even know when the Pop III epoch peaked. 

Ricotti, Gnedin & ShuU (2002) addressed this ques- 
tion using hydrodynamic cosmological simulations that 
included the feedback of both photo-dissociating and 
photo-ionizing radiation from primeval galaxies. Al- 
though their simulations were of even lower resolution 
than MBA and Y03 and parameterized star formation 
with ad hoc prescriptions, they found a cosmic star for- 
mation history quite similar to that predicted by Y03 de- 
spite small box size effects. They too could not simulate 
sufficiently large boxes to obtain converged results in the 
interesting redshift regime 20 > z > 6, but rather found 
that rare luminous objects strongly affected the average 
SFR within the simulation volume. It is fair to say that 
we have a much better idea when the first Pop III star 
formed in the universe than when the last one formed. 
Indeed, some models have been presented which predict 
Pop III stars forming as low as redshift 3 (Schneider et 
al. 2006). 

This review updates progress on direct numerical sim- 
ulations of Pop III star formation by several groups since 
the reviews of Bromm and Larson (2004), Glover (2005), 
and Ciardi and Ferrara (2005). New high resolution sim- 
ulations have been carried out by Yoshida et al. (2006), 
O'Shea & Norman (2007a) and Gao et al. (2007) which 
test the robustness of the basic picture described above 
by pushing the collapse to higher central densities, al- 
though not yet stellar densities, and by considering en- 
sembles of simulations. The effects of a FUV back- 
ground on Pop ni star formation has recently been revis- 
ited by O' Shea & Norman (2007b) and by Wise and Abel 



(2007). It is found that H2 cooling remains important in 
halos with masses approaching 10^ Mq and FUV back- 
grounds as strong as ^21 = 1- This suggests that the Pop 
III epoch may be more extended than previously thought 
and therefore occurring in rather different environments 
than originally simulated. 

Pop in star formation has been studied in new and dif- 
ferent environments as well. Pop III star formation in 
warm dark matter models has been studied by O'Shea 
& Norman (2006) and by Gao & Theuns (2007). Finally, 
Pop III star formation in gas that has been pre-processed 
by earlier generation of Pop III stars has been simulated 
by O'Shea et al. (2005), Mesinger, Bryan & Haiman 
(2006), and by Yoshida et al. (2007a,b). The surprising 
result is that despite the diversity of the environments 
studied, it is found that Pop III. 2 stars form in basically 
the same way as Pop III.l stars and that they are mas- 
sive, although they appear to be somewhat less massive 
than Pop III.l based on lower cloud core temperatures 
and accretion rates. However, the new simulations show 
that environment and formation redshift/history can have 
a substantial effect on protostellar accretions rates, sug- 
gesting some spread in the primordial IMF. 

No 3D simulation has yet been done that follows the 
entire accretion history of the star, and therefore final 
masses are still uncertain. However, approximate inte- 
grations have been done assuming spherical symmetry 
which suggest that the entire protostellar envelope can 
be accreted (Omukai & Palla 2003, Yoshida et al. 2006, 
2007a,b). Verifying this result with fully 3D simulations 
represents a grand challenge for the field. 

2. FORMATION OF POP III.1 IN ACDM: 
NEW RESULTS 

Hereafter we will refer to Pop III stars forming from 
undisturbed cosmological initial conditions Pop III.l, in 
keeping with the taxonomy introduced by McKee (these 
proceedings). Our understanding of the formation of Pop 
III.l stems primarily from high resolution 3D hydrody- 
namic cosmological simulations supplemented with pri- 
mordial gas chemistry including the all-important reac- 
tion kinetics of molecular hydrogen. The simulations of 
Abel, Bryan & Norman (2000, 2002) and Bromm, Coppi 
& Larson (2002) showed that sufficient H2 is formed in 
the cores of mini-halos of mass 10^ — lO^M© at z ~ 20 
so that ro- vibrational line cooUng of H2 will precipitate 
Jeans instability. Typical densities and temperatures in 
this gas prior to collapse are n ^ lO^cm^^ and T ~ 200 
K, implying a Jeans mass of ~ 500 M©. 

The collapse of this gas initially proceeds quasi- 
statically, and then dynamically when the density thresh- 
old for 3-body H2 formation is reached at n ~ lO^cm"^. 
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FIGURE 1. Formation of a primordial protostar using an 
ultra-hiigh resolution cosmological SPH simulation. Projection 
of gas density on various scales. The central 1 pc region at top 
right contains ~ 300 Mp , of gravitationally unstable primordial 
gas cooling by molecular hydrogen line radiation. The dark 
patch at lower left is a 0.01 M0 fully molecular hydrogen, 
optically thick cloud core cooling by coUisionally induced 
emission with central density 3 x lO'^ cm^-'. Accretion of the 
massive envelope onto this seed will form a Pop III star at 
z 19. From Yoshida et. al. (2006). 



The collapse proceeds from the "inside out", analogous 
to protostellar cloud collapse models that have been stud- 
ied for decades in connection with current day star for- 
mation as reviewed by Larson (2003). The inf ailing en- 
velope assumes a density distribution that is well approx- 
imated by a powerlaw: n[r) ^ ^^2.2 although the tem- 
perature distribution is distinctly not isothermal. Both 
groups found that the collapsing core did not fragment 
into smaller units, implying a massive star would even- 
tually be formed. Both the Abel and Bromm simula- 
tions were forced to stop well before stellar density was 
reached due to limits of physics and resolution. The sim- 
ulation by Abel, Bryan and Norman (2002) achieved a 
central density of n ~ lO'^cm^^ on scales 100 AU 
where a fully molecular core of ~ 1 M© was found ac- 
creting its envelope at a rate of ~ lO^^Mo/yr. At this 
rate >100 Mf , of material could accrete within a Kelvin 
time, suggesting that Pop III.l stars were very massive. 

New simulations have confirmed and extended this ba- 
sic picture. Table 1 summarizes the simulation details 
for the adaptive mesh refinement (AMR) simulation of 
O'Shea & Norman (2007a; hereafter ON07a) and the 
smoothed particle hydrodynamics (SPH) simulation of 
Yoshida et al. (2006; hereafter Y06). Both groups as- 
sumed WMAP 1 cosmological parameters, had similar 
mass and spatial resolution, and ran their simulations un- 
til a central density of lO'^ particles/cc was reached. The 



SPH simulation included more physics and is thus more 
reliable in the range of densities where radiation trans- 
port effects, collisionally-induced-emission (Ripamonti 
& Abel 2004), and equation of state effects become im- 
portant (n > lO'^cm-^). Y06 ran a single simulation in 
a box 0.3 Mpc on a side, while ON07a ran a dozen sim- 
ulations in three different box sizes: 0.3, 0.45 and 0.6 
/j^^Mpc. We will therefore use the ON07a simulations 
to discuss the effect of formation redshift and environ- 
ment, the Y06 simulation to discuss high density evo- 
lution, and compare the two of them over the range of 
densities where they can be compared. 



2.1. Comparison of AMR and SPH results 

Fig. 1 shows the radial density and temperature pro- 
files when the central density has reached lO'^cm^^. Al- 
though the ON07a simulation is not valid at these high 
densities, the density profiles agree extremely well over a 
large range of radii and enclosed masses considering the 
very different numerical techniques employed. The den- 
sity profile is well approximated as power law n ^ r^^-^ 
over many decades in radius. The temperature distribu- 
tions also agree extremely well for n < 10"^cm^^ where 
opacity effects become important. The temperature dis- 
tribution reflects four regions. Proceeding from larger 
to smaller radii we have first the region of cosmologi- 
cal infall, where the temperature increases to its maxi- 
mum value lOOOK) at the virialization shock located 
at r 100 pc. The second is a cooling zone where the 
gas temperature drops to the minimum allowed by H2 
cooling (r-^ 200 K) at a radius of a few pc. The third sees 
the gas temperature rise again to ^ lOOOK due to com- 
pressional heating in the quasi-statically collapsing cloud 
core at r ~ .01 pc. The fourth and innermost (r < 10^^^ 
pc) is a cooling zone accompanying the formation of 
the fully molecular core due to the 3-body H2 reaction. 
This innermost cooling zone is absent in the Y06 simula- 
tion because the cooling radiation is trapped at densities 
n > lO'^COT^^, whereas it is allowed to escape unphysi- 
cally in the ON07a simulation. When opacity effects are 
included, the central temperature rises to 2000 K. The 
difference in thermal structure seems to have little effect 
on the density profile, however 

Fig. 3 shows a comparison of profiles for molecular 
hydrogen fraction and radial velocity, analogous to Fig. 
2. Three body H2 formation creates a fully molecular 
core of size lO^'* pc and mass ~ IMq. The massive 
envelope accretes onto this seed, driving its thermal and 
chemical evolution. A maximum infall velocity of —4 
km/s is achieved in both simulations at a radius of ^ 10^^ 
pc. 

The collapse of the cloud core is relatively unimpeded 
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TABLE 1. Simulation parameters 



AMR I SPH 



Code 


Enzo 


GADGET 


Cosmological model 


WMAPl 


WMAPl 


Box size (h~^ Mpc) 


0.3, 0.45, 0.6 


0.21 


Simulations 


12 


1 


mdm(M0) 


2.6, 8.8, 21 


0.1 


mgas(M©) 


0.4 -» 5 X lO"'' 


0.015 




1015 


1015 


Physics 


H, He, H2 chemistry 
optically thin cooling 
(H2 only) 


H, He, H2, HD chemistry 
H2, HD, and CIE cooling 
radiative ttansfer 


Minimum length 
scale (pc) 


3 X 10"^ 


3 X 10"^ 


Minimum mass 
scale (M0) 


10-3 


10-2 






FIGURE 2. Radial density and temperature structure of Pop III star forming cloud, as foimd from AMR (O'Shea & Norman 
2007a) and SPH (Yoshida et al. 2006) numerical simulations. Note the excellent agreement for n < I0^^cm~^{r > 10~^ pc) where 
the gas is optically thin to H2 cooling radiation. 
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FIGURE 3. Radial profiles of H2 fraction and radial velocity from the AMR and SPH simulations referred to in Fig. 2. 



by angular momentum effects. Fig. 4 shows the spe- 
cific angular momentum profiles from both simulations. 
Again, the agreement is excellent both qualitatively and 
quantitatively. The star forming cloud has a sharply in- 
creasing specific angular momentum distribution as a 
function of enclosed mass, which is well-approximated 
by L{M) oc M. This implies the first stars form from very 
low angular momentum material. A plot of the ratio of 
the circular velocity to the local Keplerian velocity on 
mass shells (dashed Une of Fig. 4b) indicates that the in- 
falling envelope is never rotationally supported. 

How is this angular momentum profile established? 
The multiple Unes plotted in Fig. 4a indicate different 
output times corresponding to the density profiles shown 
in Fig. 2. The change with time of the angular momen- 
tum profile implies redistribution of angular momentum 
within the innermost 1000 Mq of the cloud. This is seen 
in both AMR and SPH simulations, and is understood to 
be the result of angular momentum segregation, wherein 
lower angular momentum fluid elements in the turbu- 
lent cloud preferentially sink to the center (O'Shea 2005, 
Y06). 



2.2. Evolution to higher densities: 
fragmentation? 

At densities above 10*^ cm^^ 3-body formation of H2 
becomes efficient at converting the mostly atomic gas 
into molecular hydrogen. At densities above lO'^cm^^ 
it is necessary to include the trapping of the H2 line ra- 
diation, collisionally induced emission (CIE), and chem- 
ical heating of the gas as the core becomes fully molec- 
ular. These effects have been included in the simulation 
of Y06, and therefore we refer to their results in the fol- 
lowing. 

The effect of radiation trapping and chemical heat- 
ing is to raise the temperature of the core to ^ 2000 K 
as n reaches lO^^cm"^. Prior to reaching these densi- 
ties and temperatures the core passes through a chemo- 
thermally (CT) unstable regime which Silk (1983) cal- 
culated would lead to fragmentation. Y06 evaluated 
the growth time of the CT instability using the con- 
ditions in the center of their cloud, and found that it 
briefly becomes comparable to the dynamical time at 
n ~ lO'^cOT^^ (Fig. 5), but that there is insufficient ma- 
terial satisfying these conditions to form multiple frag- 
ments. 

CIE becomes the dominant coolant above 
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FIGURE 4. Specific angular momentum profiles in Pop in 
star forming cloud. Top: time evolution in the AMR simulation. 
Bottom: final time in the SPH simulation 



n ~ lO^'^cm"^, however the dynamical timescale in 
the core remains short compared to the timescale for 
isobaric modes of the thermal instability, and therefore 
the core does not fragment into multiple objects. As 
the central density and temperature increases, the Jeans 
mass continues to drop, suggesting low mass stars could 
be produced if a fragmentation mechanism existed that 
operated quickly enough. However the simulations and 
the Unear analysis reveals that multiple fragments are 
not produced, but rather a single fragment which is the 
accretion center for the rest of the cloud. Abel refers to 
this process as "Russian doll fragmentation." This result 
needs to be kept in mind when interpreting the claims 
that very low mass stars can results from dust-induced 
fragmentation (Ciardi, these proceedings). 




10^° W 

n [cm"] 



FIGURE 5. Evolution of chemo-thermal instability in the 
collapsing Pop 111 cloud core. A growth parameter >2 is re- 
quired for the core to fragment into multiple pieces. The core 
does not fragment, but produces a single protostellar seed 
which accretes the massive envelope. See Yoshida et al. (2006) 
for more details. 



2.3. Evolution of the protostar 

3D simulations have not yet been able to evolve the 
cloud core to stellar densities, although valiant attempts 
are being made as reported by Yoshida and Turk at 
this meeting. Additional physical processes must be in- 
cluded, such as coUisional dissociation of H2 and the 
trapping of continuum radiation due to the CIE opac- 
ity. The further evolution of the core depends critically 
on the role of angular momentum in halting core col- 
lapse and building a centrifugally supported disk through 
which most of the envelope mass would presumably be 
accreted (Tan & McKee 2004). In the absence of rota- 
tion, we have the ID spherically symmetric models of 
Omukai and Nishi (1999) and Omukai & Palla (2003), 
which show that the entire envelope can be accreted pro- 
vided the late time accretion rate does not exceed the Ed- 
dington rate of m = 4 x lO^^Mo/yr. A plot of the instan- 
taneous accretion rate at the end of the Y06 simulation is 
shown in Fig. 6-top. The accretion rate drops below the 
Eddington-limited rate after 10 M© has been accreted, 
just as the protostar begins its Kelvin-Helmholtz contrac- 
tion to the ZAMS (Fig. 6-bottom). Thus most of the pro- 
tostar' s pre-main sequence evolution is sub-Eddington, 
implying that rotation, not radiation, sets its ZAMS mass. 
The three hues in Fig. 6b show protostellar evolution 
tracks for different assumptions about how much cen- 
trifugal forces reduces the accretion rate onto the proto- 
star. For even quite significant reductions in the accretion 
rate by rotation, the Pop III star is predicted to enter the 
main sequence with >60 M©. 
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FIGURE 6. Top: Instantaneous accretion rate onto the pro- 
tostar depicted in Fig. 1. Bottom: Inferred evolution of the pro- 
tostar based on a spherically-symmetric evolution model for 
1, 2/3, and 1/3 the accretion rate shown above. Dots indicate 
cessation of contraction when CNO cycle hydrogen burning 
begins. From Yoshida et al. (2006). 



2.4. Dependence on formation redshift and 
environment 

ON07a and Gao et al. (2007) have studied the de- 
pendence of the central protostellar cloud's properties 
on formation redshift and cosmological environment. As 
pointed out by White and Springel (2000), dark matter 
mini-halos in the mass range capable of forming Pop 
III.l would be expected to form in "protocluster" en- 
vironments at redshifts considerably higher than 20. A 
pure dark matter N-body study of such an environment 
by Gao et al. (2005) found that formation redshifts as 
high as 50 are likely (assuming CTg ~ 0.9). In this case, 
the higher densities and temperatures in the collapsing 
cloud would shorten formation timescales and increase 
accretion rates. ON07a simulated 12 cases by varying 
box size and random seed assuming a WMAPl cosmo- 



FIGURE 7. Density and temperature profiles for an ensem- 
ble of 12 AMR simulations when central densities have reached 
lO^'cm^-'. Density profile is insensitive for formation redshift 
and environment, while temperature profile shows considerable 
variation in the central I pc. From O'Shea & Norman (2007a). 



logical parameters. Their sample of halos formed stars 
over the redshift interval 33 > z > 19. Gao et al. (2007) 
simulated 8 cases using a multimass technique whereby 
a select region from a much larger cosmological volume 
is re-simulated at much higher mass and spatial reso- 
lution. They considered WMAPl and WMAP3 cosmo- 
logical parameters (differing principally in the amount 
of small-scale power), and also different baryon frac- 
tions. They were thus able to study halos forming stars 
over the larger redshift interval 50 > z >10. Both simula- 
tions terminated their evolutions when central densities 
of 10^'' cm^^ were reached and ignored the effects of a 
photo-dissociating UV background and opacity effects in 
the collapsing cloud core. 

Fig. 7 shows the density and temperature distributions 
at the end of the 12 simulations in ON07a. The density 
profiles are very similar to one another despite differ- 
ences in formation redshift and environment, while the 
temperature profiles show a considerable amount of vari- 
ation in the core region. This is reflected in the instanta- 
neous accretion rate plots (Fig. 8). While all 12 simula- 
tions show the generic feature of an early rapid rise to 
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FIGURE 8. Implied accretion rate versus time from an en- 
semble of 12 AMR simulations of Pop in star formation, as 
inferred from the instantaneous accretion rate at the end of each 
simulation (when central densities reached lO'^cm"^). From 
O'Shea & Norman (2007a). 




FIGURE 9. Effects of gravitino warm dark matter on the 
mass scale where perturbations are suppressed by free stream- 
ing. Cosmological mini-halos hosting the formation of Pop III 
stars become suppressed for mwdm < 20 keV. From O'Shea & 
Norman (2006). 



a maximum followed by a power-law decline, there is a 
two order of magnitude scatter in peak accretion rates. 
ON07a found a correlation between redshift of collapse 
(which correlates with box size) and peak accretion rate, 
with higher redshift minihalos exhibiting lower accretion 
rates. This was understood to be the result of the early 
rapid formation of H2 in the denser, hotter high redshift 
mini-halos which results in lower final core temperatures 
and hence accretion rates, since m ~ / G. 

Gao et al. (2007) also found considerable scatter in 
accretion rates, but no redshift correlation. Rather, they 
found their cores had a variety of shapes and amount 
of centrifugal support, and that their centrifugally- 
supported disk-like cores had the lowest accretion rates. 
Clearly, both of these effects could be operating and 
bear further study. Given the importance of the late-time 
accretion rates on setting the stellar mass, a careful study 
of angular momentum evolution in the cloud cores is 
needed. 



3. FORMATION OF POP III.1 IN AWDM 

O'Shea & Norman (2006; hereafter ON06) and Gao & 
Theuns (2007; hereafter GT07) have considered the for- 
mation of Pop III star in a universe whose dominant mass 
component is warm dark matter (WDM). WDM has been 
advanced as a solution to the many problems of CDM on 
sub-Mpc scales (e.g.. Bode, Ostriker & Turok 2001). Op- 
erationally, WDM simulations are similar to CDM simu- 
lations except that the CDM power spectrum is assumed 



to by exponentially suppressed below some scale set by 
the free streaming scale of the WDM particle, which de- 
pends upon its mass for thermal relics. Fig. 9 shows the 
suppression mass scale for gravitino WDM, given by 




(1) 

As can be seen, the mini-halos which form Pop III. 1 
stars in ACDM become suppressed for mwDM<'2Q kev, 
raising the interesting question of what is the nature of 
the first cosmological objects in AWDM ? 

ON06 carried out 8 high resolution AMR simulations 
with resolutions similar to Table 1 however varying the 
WDM particle mass over the range 12.5 < mwdm < 40 
keV. They followed the evolution until the baryonic core 
had collapsed to baryon densities of n = lO'^cm^^. Fig. 
10 shows the collapse redshift as a function of mwdm- 
ON06 found that the mwDM=^0 keV evolution is virtu- 
ally indistinguishable from the ACDM case computed 
with the same random seed, which collapsed at z=17.8. 
The effects of decreasing the WDM particle mass is to 
modify the dark matter distribution at virialization from 
spheroidal halos to elongated filaments (Bode, Ostriker 
& Turok 2001) and to reduce the core collapse redshift. 
The mwDM=l2.5 keV case produces a collapsing primor- 
dial cloud core at z=12 in the center of a cosmological 
filament (Fig. 11) which is considerably more massive 
(~ 10^ M0) than its CDM counterpart. Despite these dif- 
ferences, it is found that the high density (n > 10^ cm~^) 
evolution of the central ~ 500M© of the collapsing cloud 
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FIGURE 10. Collapse redhsift versus WDM particle mass 
for gravitationally unstable primordial cloud core. From 
O'Shea & Norman (2006). 



core is very similar to the CDM case (Fig. 12), implying 
that similarly massive stars would be produced. 

A different conclusion was reached by GT07, who 
simulated the more extreme case mwDM='i keV. In this 
case, the suppression mass scale is 10** M0. Objects 
at this mass scale virialize into long 100 kpc co- 
moving) slender filaments with no dark matter substruc- 
ture around which the baryons coalesce. Their simula- 
tion was terminated when the filament achieved a central 
density of n ^ 10^ cm^^ and temperature T ^ 200 K, but 
before it fragments. Presumably gravitational instability 
would fragment the filament into lumps of order the cen- 
tral Jeans mass which is a few hundred solar masses for 
the values above. Yet they state otherwise. They claim 
the filament will continue to collapse to higher densi- 
ties in analogy to ACDM cloud cores, and thereby frag- 
ment into lower mass protostellar seeds. They suggest 
the seeds will then accrete from the cylindrical gaseous 
envelope while they merge along the filament axis, pro- 
ducing a broad distribution of Pop 111 stellar masses. This 
scenario needs to be investigated with high resolution nu- 
merical simulations because as we have seen, whether 
the cloud fragments or not depends not only on the lo- 
cal Jeans mass, but on a competition of timescales which 
can only be determined from a self-consistent multidi- 
mensional calculation. 



4, FORMATION OF POP IIL2: 
VARIATIONS ON A THEME 

In this section we return to simulations of the formation 
of Pop III stars in a ACDM universe which examine the 




FIGURE 11. Formation of Pop III star in a cosmological 
filament in a AWDM universe. From O'Shea & Norman (2006) 




FIGURE 12. Evolution of density profile in the collapsing 
cloud core in the mwDM=^'^-5 kev simulation shown in Fig. 11. 
For ;i > 10^ cm the ~ 500 Mp, core evolves very similarly 
to the ACDM case despite the delay in achieving collapse 
conditions. From O'Shea & Norman (2006). 
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influence of radiation backgrounds on their formation 
history. We discuss two cases in particular, the forma- 
tion in the presence of a concurrent photo-dissociating 
Lyman- Werner background, and formation from gas that 
has previously been ionized by stellar sources (relic HIT 
regions). 

4.1. Effect of a Lyman- Werner UV 
background 

Machacek, Bryan & Abel (2001) showed using mod- 
erate resolution AMR survey simulations that a Lyman- 
Werner background radiation field delays but does not 
extinguish Pop III star formation by raising the minimum 
halo mass capable of cooling by H2 line transitions. They 
carried out simulations for fixed LW backgrounds corre- 
sponding to 721 = 0, 10^^-^ 10"'^. The formula for a fit 
to their minimum halo mass is given by Ciardi in these 
proceedings, but generally rises from M,„,„ = 10^ Mg at 
721=0 to M,nin = lO'^Ma at 721 = 10"' '. Wise and Abel 
(2005) extrapolated the fit to LW backgrounds as high as 
721=1 and used a Press-Schechter-type analysis to esti- 
mate the evolution of the LW background for three values 
of the assumed mass of the the Pop 111 stars produced. 
They found that Pop III star formation becomes self 
Umiting by z=25 and that the LW background reaches 
721 = 0. 1 - 1 in the redshift range 20 > z > 10. Y03 car- 
ried out SPH survey simulations similar in resolution to 
those of MBA and found that H2 cooling becomes in- 
efficient in primordial mini-halos for 72i>0.01, making 
the extrapolation of the MBA minimum halo mass fitting 
formula to 72i=l somewhat suspect. Neither the MBA 
nor the Y03 simulations had the mass or spatial reso- 
lution to simulate baryonic core collapse. Motivated by 
this, O'Shea & Norman (2007b) carried out such simu- 
lations. 

Fig. 13 shows the time evolution of central values for 
gas density, temperature, entropy, and H2 fraction for 
four values of 72i . The figures show the time delay effect 
very clearly. The physics is very simple. For 72i>0.01 the 
photodissocation time is shorter than the 2-body forma- 
tion time, consequently the H2 fraction assumes its equi- 
Ubrium value. The cooUng time is then proportional to 
721, ™d reaches 10^ yr for 72i = l. The surprising result is 
that collapse eventually occurs even in the extreme case 
721=1 despite the greatly reduced H2 fraction in the core. 

The explanation for this result is that the time delay 
allows the virial mass and temperature of the halo to 
increase through mergers to Mv,> ~ lO^M© and T^ir ~ 
2000 K where the per molecule cooUng rate in the core 
is greatly elevated. In the most extreme case, central 
temperatures reach 5000 K which accelerates production 
and cooling of H2 to the point of baryonic collapse. 




10« 1.5xl0« 2xl0« 10* 1.6xl0« 2xl0« 

t [Myr] t [Myr] 




t [Myr] t [Myr] 

FIGURE 13. Evolution of central density, temperature, en- 
tropy, and H2 fraction in a mini-halo for different Lyman- 
Wemer UV background intensities (solid: Tzi = 10~^, short 
dash: 72i = 10~^, long dash: 72i = 10~', dot dash: J21 = 1). 
In the strongest UVB, gravitational collapse is delayed by 100 
Myr, but not prevented. From O'Shea & Norman (2007b). 

Fig. 14-top shows the increase of the mass-weighted 
temperature in the central lOOM© of gas at the collapse 
redshift for increasing values of 72i . This translates into 
larger accretion rates onto the protostar at early times 
(Fig. 14-bottom) exceeding lO^^Mo/yr for 72i>0.01. 
Taken at face value, this suggests that Pop III. 2 stars 
forming in high LW backgrounds may be more massive 
than Pop III. 1 stars, although we must bear in mind that it 
is the late-time accretion rate that governs the final mass 
of the star. This suggestion bears further study. 

Wise & Abel (2007) independently carried out a study 
similar to O'Shea & Norman (2007b), but considered 
even more massive halos (Myir < 4x10^ M©) where H2 
cooling gives way the Lyman alpha cooling. They find 
a smooth transition in collapse redshift and halo mass 
between those that cool by H2 fines and those that cool 
by atomic hydrogen lines (Fig. 15). They find that col- 
lisional ionization of H in the warmest halos produces 
enough electrons to catalyze the formation of H2 even 
in the presence of strong LW backgrounds. The conclu- 
sion from both of these works is that H2 cooUng remains 
the dominant mediator of early star formation even for 
strong UVBs, and that by restricting ones attention to 
only those halos for which r^y > lO'' K one misses a 
large fraction of star forming halos. 
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FIGURE 14. Conditions in the central 100 Mq cloud core as 
a function of LW mean intensity. Top: mean temperature, Bot- 
tom: accretion rate. Because a LW background delays baryonic 
collapse, virial temperatures, core temperatures and hence ac- 
cretion rates are elevated in late-forming primordial mini-halos. 
From O'Shea & Norman (2007b). 
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FIGURE 15. Smooth transition from H2 coolers to Lyman 
alpha coolers, in terms of collapse redshift and central temper- 
ature. From Wise & Abel (2007). 



4.2. Formation in a relic HII region 

Finally we consider the case of a mini-halo that has 
previously been ionized by a strong UV source, Pop III 



or otherwise. As reviewed by Ciardi (these proceedings). 
Pop III stars photo-evaporate most of the baryons from 
the halos that formed them, and create large HII regions 
of 1-2 kpc in diameter. After a few Myr the Pop III star 
dies, leaving behind a relic HII region. Gas cools and 
recombines out of equihbrium, and the large residual 
electron fraction quickly catalyzes the formation of H2 
via the H- reaction (O'Shea et al. 2005). The question 
is under what conditions does a neighbor halo form 
a primordial star, and in those that do, what are their 
properties? 

The survival of neighbor halos depends on a number 
of factors including halo mass, evolutionary stage, prox- 
imity to the ionizing source, and its ionizing luminosity 
(Ahn & Shapiro 2007, Whalen et al. 2007). The most 
important factor appears to be the central density of the 
neighboring halo at the time of the time it is illuminated 
by the ionizing source. Densities of n ~ 10^ cm~^, corre- 
sponding to the "loitering phase" of the evolution of pri- 
mordial cloud cores, are high enough to self-shield the 
cloud against the LW radiation from the first star. Cores 
of lower density are photo-dissociated, but crashed to 
higher density by the cometary 1-front. The residual ion- 
ization in these cores rapidly produces H2 after the first 
star dies. Only the lowest density neighbor halos have 
their baryons stripped, and are therefore unable to form 
stars. Overall, the effect on star formation appears to be 
neutral (Ahn & Shapiro 2007). 

O'Shea et al. (2005) carried out a high resolution 3D 
AMR cosmological simulation of the evolution of a halo 
forming inside the relic HII region of a Pop 111.1 star. 
Resolutions similar to that given in Table 1 were used 
to follow the collapse of the second halo, which be- 
gan 23 Myr after the first star died. The collapse pro- 
ceeded in a way very similar to that of a Pop III.l star 
once central densities of IC^cm^^ were reached, indi- 
cating the universality of how Pop III stars are formed. 
However, the instantaneous accretion rate when cen- 
tral densities of 10 '"cm"^ were reached was quite low 
(m ~ lO~'*M0/yr) due to the fact that the second halo 
had a large amount of angular momentum and produced 
a centrifugally-supported disk in its center. 



4.3. Effect of preionizadon and HD cooling 

So far we have considered only Pop III stars forming 
from gas that has remained neutral since recombination. 
Stars forming in chemically pristine gas that has previ- 
ously been ionized start from initial conditions that dif- 
fer thermally, chemically, stracturally, and kinematically 
from Pop III.l. Photoionization will raise the entropy of 
the gas and smooth out the clumping due to photoheating 
and photo-evaporation (Oh & Haiman 2003). In addition. 
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the high electron fraction will catalyze rapid formation of 
H2 molecules via the H- channel when the HII region be- 
gins to recombine after the ionizing source switches off 
(O'Shea et al. 2005). In the absence of a high LW back- 
ground, H2 fractions are elevated by 2 orders of magni- 
tude at the cosmic mean density relative to the Pop III.l 
case. Finally, these conditions are ideal for the formation 
of HD, which becomes an important coolant at moder- 
ate densities (10 < n < lO^cm^^) and low temperatures 
(T<150K) due to chemical fractionation and HD's large 
permanent dipole moment (Glover, these proceedings). 

The effect of these processes on the evolution of Pop 
III star-forming halos has been examined by Mesinger, 
Bryan & Haiman (2006) (see also Bryan et al., these pro- 
ceedings). They carried out high resolution AMR sim- 
ulations using Enzo of a region that forms a Pop III.l 
star at z=25 in the absence of radiative feedback. With- 
out HD cooling, they found that the negative and positive 
feedback effects roughly cancel, with negative feedback 
(density smoothing) slightly outweighing positive feed- 
back (accelerated H2 formation). Pop III. 2 stars formed 
in such regions are found to have very similar properties 
as Pop III.l stars as estimated by mass accretion rates. 
When HD cooling is included, however, the simulations 
predict a reduced mass accretion rate because of the 
lower tempersture and hence soundspeed in the accreting 
gas. Pop III.2 stars with typical masses of 30 — 40 M© are 
claimed. 



4.4. Second generation primordial stars 

A similar conclusion was reached by Yoshida et 
al. (2007a,b), who recently have simulated the second 
primordial star to form in a given primordial mini- 
halo. They self-consistently simulate the formation of a 
Pop III.l star as in Y06, and follow the expansion of 
the HII/Helll region it produces. As in earlier studies 
(Whalen et al. 2004, Kitiyama et al. 2004), gas in the 
minihalo core is swept out to the virial radius by a strong 
photo-evaporative flow. After the Pop III.l star dies, this 
gas recombines and cools, forming HD. They continue 
the calculation until this gas recollects into the dark mat- 
ter mini-halo core and forms a Pop III. 2 star. They find 
that it takes about 100 Myr for the second primordial star 
to form, and that as a consequence of HD cooling, the 
mass accretion rate is reduced, in agreement with the re- 
sults of Mesinger, Bryan & Haiman (2006). In Yoshida, 
Omukai & Hemquist (2007) they calculate the accretion 
history of the Pop III. 2 protostar using a ID protostellar 
accretion code. The result is shown in Fig. 16. It enters 
the main sequence with a mass of ~ 30 M©, significantly 
less than the Pop III.l star that pre-processed its birth 
cloud. 




0.1 1.0 10.0 100.0 
protostellar mass [Mg] 



FIGURE 16. Evolution of the protostellar mass and radius 
for a Pop III.2 star forming in the halo of an earlier Pop III. 1 star 
which photoionized its environment. From Yoshida, Omukai & 
Hemquist (2007). 
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